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By using the supersymmetry method we derive an explicit 
expression for the parametric correlation function of densities 
of eigenphases 9 a of the S-matrix in a chaotic quantum sys- 
tem with broken time-reversal symmetry coupled to continua 
via M equivalent open channels; a = 1,..,M .We use it to 
find the distribution of derivatives of these eigenphases over 
the energy ("phaseshift times") as well as over an arbitrary 
external parameter. We also find the parametric correlations 
of Wigner-Smith delay times. 

PACS numbers: 05.45. +b, 24.30 v 



The most fundamental object characterizing the pro- 
cess of quantum scattering is the unitary S-matrix relat- 
ing the amplitudes of waves incoming onto the system 
and the amplitudes of scattered (outgoing) waves. 

Recently there was a considerable growth of interest in 
the phenomenon of quantum chaotic scattering ( see re- 
), both theoretically 0^^R|,and experimen- 
Due to the chaotic nature of the underlying 
scattering dynamics the S-matrix characteristics behave 
in an irregular way when parameters of either incoming 
waves (e.g. the energy) or characteristics of the scatter- 
ing region (e.g. the form and strength of the scattering 
potential, the strength of the magnetic field through the 
ballistic microstructure, etc.) are slightly changed. Be- 
cause of this fact it seems to be most adequate to describe 
such a behaviour in terms of some statistical measures: 
distributions and correlation functions. 

At present, there are two complementary theoretical 
tools employed to study open quantum systems whose 
closed classical counterparts demonstrate chaotic be- 
haviour. These are the semiclassical approach and 
the stochastic approach Q-^]. The latter one is based 
on the notion of the universality of the chaotic scattering 
phenomenon. 

It is well known that the majority of "closed" chaotic 
quantum systems of quite different microscopic nature 
shows a great degree of universality on the level of its 
statistical characteristics of spectra and eigenfunctions 
stemming from the very basic fact of chaotic internal 
motion. Because of this universality one achieves the 
correct description of the properties of such systems by 
exploiting the similarity with ensembles of large random 
matrices (RM). The latter principle is a commonly ac- 
cepted one in the domain of Quantum Chaos ]10| , see 
also [ [ill , and recently some sound arguments were given 



in favour of its general validity |T^] . 

Provided the statistics of the system Hamiltonian is 
specified, one can work out the S-matrix by standard 
methods in the theory of quantum scattering Jl3| and 
study its statistical characteristics. This way was pi- 
oneered by Verbaarschot et al. || who calculated the 
correlation function of S-matrix elements. Other char- 
acteristics of the S-matrix and related quantities can be 
efficiently studied in the framework of this approach as 
well 

One can also try to make use of the expected universal- 
ity directly on the level of S-matrix without any reference 
to the system Hamiltonian. Such a method was devel- 
oped in a great detail in a series of papers by Mello and 
collaborators [0]. The probability density for the whole 
S-matrix can be obtained if one makes the assumption 
of minimal information content of such a distribution re- 
specting the requirements of S-matrix unitarity and sym- 
metry constraints. Provided all the relevant information 
about the system is encoded in a value of the average S- 
matrix (S) the probability P(S) was shown to be given by 
a so-called Poisson's kernel [||. More recently, Brouwer 
p7[ has shown that the general Poisson kernel can be de- 
rived directly from the Hamiltonian approach if one takes 
the Hamiltonian matrix from a quite specific Lorentzian 
RM ensemble. In view of the asymptotic equivalence of 
the Lorentzian and Gaussian matrices |T^ ] it is natural to 
expect that the Poisson's kernel is of the same universal- 
ity as the spectral statistics for closed systems discussed 
above. 

If one wishes to study the dependence of S-matrix on 
external parameters without explicitly considering the 
system Hamiltonian, one should make some additional 
statistical assumptions beyond the minimum information 
approach. One possible way is to simulate such a depen- 
dence by a kind of " Brownian motion" in the correspond- 
ing S-matrix space. Quite a detailed presentation of such 
a method with some attempts to find a microscopic jus- 
tification of the whole procedure can be found in |L8[ . 
It turns out, however, that the Brownian motion pic- 
ture is in disagreement with the results obtained starting 
from the Hamiltonian formalism. Therefore, the latter 
approach seems to be the only consistent one when we 
are interested in the parametric variations. 

Apart from the S-matrix elements, the scattering phase 
shifts 6 a (defined via the S-matrix eigenvalues expi9 a ) 
are used intensively to characterize the chaotic scattering, 
see e.g M. Quite recently their statistical characteris- 
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tics were studied numerically in some detail |1S|,|20[| . The 
derivatives of phase shifts over the energy r a = 80 a /dE 
( which we suggest to call the "phaseshift times" ) are 
particularly interesting and related to the mean time 
spent by a particle in the interaction domain (Wigner- 
Smith delay time, see p2| and references therein) as t w — 
jj^T a . The fluctuations of the individual phaseshift 
times r a can be therefore used to characterize the possible 
fluctuations of delay times in the scattering process. Let 
us mention that various aspects of the general problem of 
delay time characteristics is under quite intensive study 
for a while JT5| , pl|j2^ ] . More general parametric deriva- 
tives of the scattering phase shifts can also be related to 
some observable quantities. As a particular example we 
mention the relation between the persistent currents and 
the derivative of the total phase shift over the magnetic 
flux ||. 

In order to be able to study statistics of phase shifts 
9 a within the Hamiltonian approach, it is convenient to 
use the following representation for the scattering matrix 
(see e.g. @): 

S= [i-iirW*R(E)W] x [I + iirW^ R(E)W}~ 1 (1) 

where we introduced the resolvent R(E) = (E—H)^ 1 cor- 
responding to the Hamiltonian TL of the " closed" chaotic 
system and W is the matrix that couples the chaotic re- 
gion with the incoming/outgoing waves. In the stochastic 
approach the system Hamiltonian TL is replaced by the 
N x N matrix Hij taken from the corresponding Gaus- 
sian RM ensemble. In order to allow for the dependence 
of the S-matrix on the external parameter X to be taken 
into account it is convenient to consider the set of RM 
H i:j to be of the form Q : H(X) = H + ^=#1- To 
be specific, we consider in the present paper the simplest 
possible case when Hq runs over the Gaussian Unitary 
Ensemble (GUE) and Hi is an arbitrary but fixed matrix 
from the same ensemble. Physically this case is known 
to describe systems with completely broken time reversal 
symmetry (TRS). The scattering in systems with partly 
or completely broken TRS is under quite intensive theo- 
retical lIQlllEil and experimental fH study. 

The matrix W is N x M matrix of amplitudes Wi a , a 
1 . 2, M, coupling the internal motion to M open chan- 
nels. Without much loss of generality these amplitudes 
can be chosen fixed in a way ensuring that the average 
S— matrix is diagonal in the channel basis: S a b = S a bS aa 
H . Provided the energy E is real, one finds the following 
expression p|: 



~n 1 " 19(E) 



Saa 



1 + 19 {E) ' 



Wi, 



(2) 



where g(E) = iE/2 + (1 — E 2 /4) 1 / 2 and we assumed 
that all M channels are statistically equivalent for the 
sake of simplicity. The strength of coupling to continua 
is convenient to be characterized via the " sticking proba- 
bilities" (also called the " transmission coefficients" [f|,U ) 



T a = 1 — \S aa \ 2 which are given for the present case by 
the following expression: 



2Re g{E) 



(7 + 7" 1 ) 



(3) 



The quantity T a measures the part of the flux in channel a 
that spends substantial part of the time in the interaction 
region . We see that both limits 7 — > and 7 — > 00 
equally correspond to the weak effective coupling regime 
T a -C 1 whereas the strongest coupling (at fixed energy E 
) corresponds to the value 7 = 1. The maximal possible 
coupling corresponding to the upper bound T a = 1 is 
achieved in the present model for an energy interval in 
the vicinity of the center E = 0. 

The eq.(jl|) shows that eigenphases 9 a considered mod- 
ulo 2-7T (2(| are determined in a unique way by the eigen- 
values z a (E, X) of the matrix A X (E) = nW*R x (E)W in 
view of the relation: 9 a = —2 arctanz a . To this end let us 
introduce the density Pe,x( z ) — jr* Sa=i ^( z ~ z a(E, X)) 
and consider the correlation function: 

{PEfi{ z l)PE+n,x{z2)) ~ {pEfl{ z l))(pE+i\x{z2)) 

where the angular brackets stand for the averaging over 
the RM ensemble. This correlation function can be easily 
found provided the following one is known for e — ► + : 



f(zi,z 2 ) = ( Tr- 
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A somewhat unpleasant feature of eq.(J5|) is that the 
random matrix H enters it only via the matrix A. How- 
ever, due to the identity 



det(z - A X (E)) = 
det[R x (E)(E - H(X) - 



(6) 



the right hand side of eq. (||) can be written in the follow- 
ing form: 
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where we introduced the notations: Z 



(p) 



J D 



P 



1,2 and H eff {X-zf) 



H(X) + 



-i(-l) p e+ 



In the latter form the expression is very close in its struc- 
ture to that used to study parametric correlations in en- 
sembles of large RM |yj and can be evaluated by es- 
sentially the same method. The general strategy is as 
follows: 1) to represent determinants in the denominator 
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of the preceeding expression by auxilliary Gaussian inte- 
grals over commuting variables and those in the numera- 
tor by similar integrals over anticommuting variables. 2) 
To perform the averaging over the GUE distribution of 
Ho. 3) After employing the Hubbard-Stratonovich trans- 
formation to exploit the limit N ^> M when integrating 
out some (" massive" ) degrees of freedom in the saddle- 
point approximation. After this is done the integral over 
the remaining ("massless") degrees of freedom can be 
represented in a form of the so-called zero-dimensional 
supermatrix nonlinear er— model introduced for the first 
time by Efetov [§7| . 

For the connected part of the correlation function 
i), we find after the set of standard manipulations 
the following expression: 



f( z i,z 2 ) 
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where we put for simplicity the energy parameter E = 
and made the natural rescaling [[fl|: ui = nil/ A; x 2 = 
(nX/NA^jjTrH 2 , with A oc l/N being the mean level 
spacing of the GUE at E = 0. 

The eq.(Q) describes the parametric correlations of 
eigenvalues of the matrix Ax (E) in closed form and pro- 
vides the basis for extracting the statistical properties of 
scattering phase shifts. First of all, it is evident that the 
integrations in eq.(Q) can be performed explicitly when 
x = lo = thus giving the expression for the correlation 
function of phase shifts at fixed value of the energy E 
and the parameter A. In order to represent it in a sim- 
ple form for any value of the coupling constant 7 it is 
convenient to introduce new variables 4> related to the 
scattering phase shifts 9 as 4> — arctan{7 _1 tan (9/2)}. 
A direct calculation shows that the correlation function 
of the densities of the angles <fi is equal to 



IC((f>l, (f> 2 ) 1^02 



sin M((f>i — (f> 2 
irM sin ((hi — d> 
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for any number of open channels. Remembering the re- 
lation between 7 and (S), see eq.(||), one can satisfy one- 
self that eq. (|J) exactly coincides with the pair correlation 
function following from the Poisson's kernel distribution 
as defined in |BiL7||. 

Let us turn our attention to the distributions of the in- 
dividual phaseshift times r a = 89 a /dE and general para- 
metric derivatives of phase shifts. In view of the relation 
9 a = — 2arctanz a these distributions can be found eas- 
ily provided the joint probability densities VE(z,t) and 
Vx{z,y) of z a and one of its parametric derivatives - 



t a = dz a ldE or y a = dz a /dX - is known. The func- 
tions Ve and Vx can be determined from the correlation 
function eq.(||) as follows (see e.g. p8|): 

-Pe(z, t) = jj (E^! S(z - z a )S(t - dzJdE)) (9) 
= Mlimn^o ^C E ,n,x=o{ z i = z, z 2 = z + tCl) 

Vx(z,y) = jf (Ea=i s ( z - z a)Kv - dz a /dX)^ , 1Q , 
= Mlimx^o XK, E ,n=a,x{zi = z,z 2 = z + yX) 

Using the procedure described above we find the fol- 
lowing distribution of inverse (scaled) phaseshift times 
u = 2tt/(tA): 



r u (u)^j I (j2^[u-2n/(T a A)]j 



(11) 
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where Iq(z) stands for the modified Bessel function, 
g = 2T- 1 - 1 and T is the "sticking probability" de- 
fined earlier. Similarly, the scaled parametric derivatives 
v a = d9 a /dx are found to be distributed according to the 
probability density 



V v (v) = 



du 



(27T) 1 /: 
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where V u (u) is defined in the preceeding equation. Both 
distributions depend on the coupling constant 7 only via 
the sticking probability T rather than via the average S- 
matrix. They assume the simplest form for the " critical" 
coupling T = 1 corresponding to the most strong over- 
lap of individual resonances allowed for the few-channel 
scattering |lj|l(| . Under this condition one finds the fol- 
lowing distribution of scaled phaseshift times t s = 1/u: 
VAT a ) = a=&Tr M -*e-VT: 

In the recent paper pCj this formula was found to be in 
a good agreement with results of direct numerical simula- 
tions. The authors also found a way to arrive at such an 
expression exploiting some plausible assumptions about 
the statistics of S-matrix eigenphases and their deriva- 
tives. 

Having at our disposal the exact distribution eq. ( |TT| ) it 
is instructive to calculate the mean value and the variance 
of the phaseshift times. One finds: 



(r) 



2tt ^ 
M~A' 



(t 2 ) - (r) 2 2M(T- 1 - 1) + 1 



M - 1 



(13) 



The first of these relations is quite well known |2l|,|22j . 
It shows that the mean delay time (r) is determined by 
the mean level spacing A of the closed system and the 
number M of open channels. On the other hand the 
magnitude of delay time fluctuations measured by the 
relative variance of the phaseshift time distribution, see 
eq.(|13), is determined both by M and T. Generically, 
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the fluctuations are the weaker the larger is the number 
of open channels M and the stronger is the coupling to 
continua: 1 — T <C T. Let us also mention as an inter- 
esting feature the divergency of the time-delay variance 
at M = 1, which is a consequence of the power law tail 



--M-2 



typical for the distribution "P t (t s ). 



Eq.(^) can be used to study the parametric correla- 
tion functions of the Wigner-Smith delay times t w (E, X) . 
However, this calculation turns out to be quite lengthy 
and will be presented in a more extended publica- 
tion. On the other hand, remembering the definition 
tw = — jj'gjs ]ndetS(E) one finds the relation t w — 
-jjlm T^E - H + inWW^)- 1 which follows directly 
from eq.(Q) upon using the identity eq.((|). Now it is 
evident that the calculation of the correlation of fluctu- 
ations of Wigner-Smith delay times Stw(E,X) = tw — 
(tw) amounts to evaluating the average product of the 
resolvents of the non-Hermitian effective Hamiltonians 
Ti ± iirWW + . For the case of chaotic systems with pre- 
served TRS such a calculation was done earlier in p2| . 
For the present case we find: 



(5t w (E,0)Stw{E + Q,X)) = 
x dX d\i cosw(Ai — A) 

2 



\maJ 

exp{^(A?-A 2 )}x i±** 

2 [1 + Xig 1 

The variance of the t w is given by: 

(r w ) ~ (rw) 2 _. 2 



(14) 
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T 2 (M 2 



-[1-(1-T) M+1 ] (15) 



showing the same qualitative features ( divergencies at 
M = 1 or T — > 0) as those following from the eq.(|l3|). 

In conclusion, we studied analytically the parametric 
correlations of scattering phase shifts and the distribu- 
tion of phase shift derivatives for a chaotic quantum sys- 
tem with broken TRS. The method of calculation auto- 
matically provides the universality of the obtained results 
for a very broad class of systems which can be described 
by the mapping on the same supermatrix nonlinear a— 
model [0 . Another important feature of the method is 
that it provides the regular basis for studying a general 
case of disordered systems where the effects of localiza- 
tion modify essentially the phaseshift statistics |l9| . This 
issue as well as the extension of our results to the case of 
systems with preserved TRS are left for future research. 
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